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• We implemented a two- and three-dimensional finite difference/front trac mg 
technique to solve liquid-vapor phase change problems. The mathematical and e 
numerical features of the method were explained in great detail in our previous 
reports. Briefly, we used a single formula representation which incorporated jump 
conditions into the governing equations. The interfacial terms were distributed 
as singular terms using delta functions so that the governing equations would be 
the same as conventional conservation equations away from the interface and in 
the vicinity of the interface they would provide correct jump conditions. We used 
a fixed staggered grid to discretize these equations and an unstructured grid o 
explicitly track the front. While in two dimensions the front was simply a connec- 
tion of small line segments, in three dimensions it was represented by a connec ion 
of small triangular elements. The equations were written in conservative forms 
and during the course of computations we used regnding to control the size of 
the elements of the unstructured grid. Moreover, we implemented a coalescence 
algorithm in two dimensions which allowed the merging of different fronts or two 
segments of the same front when they were sufficiently close. We used our co e 
to S study thermo capillary migration of bubbles, burst of bubbles at a free surface, 
buoyancy-driven interactions of bubbles, evaporation of drops rapid evaporation 
of an interface, planar solidification of an undercooled melt, dendritic solidification, 
and a host of other problems cited in the reference. 

• The code was tested against different benchmark problems with known qualita- 
tive behavior or analytical solutions. For example, we solved the one dimensiona 
film boiling of a vapor layer beneath a liquid layer using our three-dimensiona 
code Initially, the liquid and vapor were at T sat (be., saturation temperature) and 
suddenly the lower wall was set to T w (be., wall temperature) which was higher 
than the saturation temperature. As a result, the vapor film started to grow an 
interface moved upward. This problem has an analytical solution (See for example 
Alexiades and Solomon, 1993) and it is good for benchmarking. We also compared 
the results of our numerical simulations with a more challenging problem given y 
Welch and Wilson (2000). Here, the liquid existed in metastable superheated con- 
dition and vapor was at saturation temperature. In contrast to the first case where 
thermal boundary layer originated from the wall and grew over the time here, t e 
thermal boundary layer originated at the interface and remained thin during e 
simulation. The thin thermal layer poses a great deal of difficulty for the numenca 
computations. In both cases, we were able to get satisfactory resu ts wi mo es 
grid resolution (i.e. 1 - 2 1 x 3 box). In another case, we compared the results of 
our thermocapillary studies with the available results in the literature (i.e. Young 
et al. 1959, Chen and Lee 1992, and Nas 1995) and found satisfactory agreement 
We did other validation tests for solidification processes and other problems and 


the detail can be found in our previous reports. 


• We modified our code in several occasions either to make it more robust or to 
introduce new features. In the original code, we used an iterative scheme based on 
Newton’s method to find the interface velocity. This algorithms was time consum- 
ing and later on we replaced it with a non-iterative algorithm which took advantage 
of some mathematical manipulations to bypass the iterations. In three dimensions, 
we first started with a point-based algorithm to operate on the front. This algo- 
rithm relied on book keeping between the front points and computations of some 
of the necessary terms such as curvature was difficult and in some cases inaccu- 
rate. We later on adopted an element-based technique for the three-dimensional 
unstructured grid which resulted in a more accurate and robust code. 

• We computed various terms for better quantification of the key parameters. 
Among them was computations of deformation factor in two and three dimensions 
and implementations of a “cluster multiple labeling algorithm” to quantify the 
formation of clusters of bubbles in large scale simulations of thermocapillary flows. 
Again, please refer to our past reports for further detail. 

• The mathematical and numerical modeling was modified so that the code can 
handle interactions of bubbles in unbounded domains. This modification enabled 
us to follow the long time evolution of bubbles in order to characterize the time- 
averaged features of the flow. In the original code, this was not possible since the 
bubbles were confined by the container walls. 

• We modified the three-dimensional code to study open surface problems such 
as film boiling and planar solidification. The code was initially tested against 
Rayleigh- Taylor instability in a flow without phase change and produced satisfac- 
tory results. Later on, (right after grant finished), it was tested against film boiling 
on a horizontal surface and produced some of the known experimental features such 
as formation of mushroom in low Jacob numbers and variations of Nusselt number 
with time in the course of film boiling. 

• We did extensive literature survey and numerous computations to gain insight 
into the problem and improve our methodology. We have already presented and 
published parts of these studies in various international conferences and journals. 
We are in the process of finishing up a few articles which we could not complete 
due to lack of computational resources. We are currently using computational 
time awarded to us through the university of Michigan to finish some of our two- 
dimensional simulations. The three-dimensional studies are made possible by an 


award from NPACI. The following list shows are past publications, pending ones 
and the one that will be submitted in the near future. 
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